Rigorous results for tight-binding networks: particle trapping and scattering 
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We investigate the particle trapping and scattering properties in a tight-binding network which 
consists of several subgraphs. The particle trapping condition is proved under which particles can be 
trapped in a subgraph without leaking. Based on exact solutions for the configuration of a 7r-shaped 
lattice, it is argued that the bound states in a specified subgraph are of two types, resonant and 
evanescent. We also link the trapping rigorous result to the scattering problem. The scattering 
features of the 7r-shaped lattice is investigated in the framework of the Bethe Ansatz. 
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I. INTRODUCTION 



Trapping and scattering of a particle is an important 
feature in many quantum information processing sys- 
tems. Due to the development of technology, the imple- 
mentation of quantum information processing in quan- 
tum systems with periodic potential, such as optical lat- 
tices [U, arrays of quantum dots 0, photonic crystal 
Q and coupled-resonator optical waveguide [H, has at- 
tracted intensive investigations. The design of quantum 
device based on these promising technologies relies on the 
particle trapping and scattering properties in a discrete 
system. A heuristic example shows that the quantum 
confinement in a discrete system is distinct from its coun- 
terpart in continuum media [5], due to the Wannier-Stark 
localization @. 

This paper focuses on noninteracting particles on dis- 
crete lattice, which is treated by tight-binding approxi- 
mation. Intuitively, the particle trapping is implemented 
by sufficient strong on-site potential as the continuous 
system. In contrast to continuum, however, different dy- 
namical properties emerge in the lattice system due to 
its distinct dispersion relation: a local wave packet can 
be confined by linear potential distribution 0] and the 
degree of spreading of a propagating wave packet can 
be controlled by judicious choice of the particle energy 
@; HI Q • Recent studies show that Fano resonance may 
be employed to construct the perfect mirror or trans- 
parency so as to control particles in a region of the lattice 
[13, Ell) Q3 y i a engineered configurations. Because of the 
numerous varieties of the possible geometry of the quan- 
tum network, we believe it is beneficial to have lattice- 
based rigorous results and exact solutions for the devise 
of a quantum device. In this paper, we show rigorously 
that the perfect particle trapping without any leakage can 
be achieved in simple tight-binding networks. This pro- 
vides a method to devise the quantum network to confine 
particles with required mode. We also link the trapping 
rigorous result to the scattering problem. This general 
finding is illustrated by a practical network consisting of 




FIG. 1: (Color online) An arbitrary graph of tight-binding 
network within part of which particles may be confined with- 
out any leakage. The graph can be decomposed into two sub- 
graphs I and I' which are connected via the coupling between 
the joint sites (i,j,h) and (i',j',h') (a). The perfect bound 
states can be formed in subgraph I when the eigen functions 
of Hi have wave nodes on all the joint sites (i,j,h), which 
are denoted by empty circles (b). The existence of more wave 
nodes (like site g) may allow multiple bound states. 



a waveguide with an embedded 7r-shaped subgraph. Ex- 
act solutions for such types of configurations are obtained 
to demonstrate and supplement the rigorous results. 
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II. RIGOROUS RESULT FOR PARTICLE 
TRAPPING 



A general tight-binding network is constructed topo- 
logically by the sites and the various connections be- 
tween them, and is also represented as a vertex-edge 
graph. Cutting off some of the connections, a graph is 
decomposed into several subgraphs. So when a particle is 
strictly trapped within a certain region of a network, one 
can say that it is confined in a specified subgraph. The 
main aim of this paper is to answer the questions of what 
kind of subgraph can trap a particle as bound state and of 
how such a subgraph scatters a particle when it is embed- 
ded in a waveguide. The Hamiltonian of a tight-binding 
network, or a graph which consists of hq subgraphs reads 
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Hi = - Y,( K fj4,i a i,j + H - c -) + 2 vf a h a i,i> 

ffim = -^{K^aj^nj +K.C.), 
i,3 

where label I denotes the Ith subgraph of Ni site, which 
subgraph is defined by the distribution of the hopping in- 
tegrals {Kfj} and on-site potentials }> and a\ - is the 
boson or fermion creation operator at the jth site in the 
ith subgraph. Here, Hi and Hi m represent the Hamiltoni- 
ans of the subgraphs and the couplings between them. In 
terms Hi m , site i (j) is the joint site of subgraph I (m) for 
the connections to other subgraphs. Obviously, the de- 
composition of subgraphs is arbitrary, and can be imple- 
mented at will. Figure [1] shows an example schematically. 
Note that the Hamiltonians Hi (also H) are quadratic in 
particle operators and can be diagonalized through the 
linear transformation 

ni, (2) 

3 

which leads to 

Hi = ^2ei,kvlkVl,k, (3) 

k 

where Si t k is the corresponding eigenvalue of Hi for the 
eigenfunction g l k ^ . Site j is defined as the wave node 

for the eigen mode k of graph I if we have g k j = 0. 
We denote the wave node as j(l,k), which reflects the 
property of the eigen state rf lk |0) of Hi 

aijvU 1°) = °> ( 4 ) 

where |0) is the vacuum state. Now we consider the case 
of that all the joint sites of the subgraph I are the wave 
nodes of eigen mode k. Under this condition, we have 

H (vL |0>) = Hi (ril ik |0)) = e Lk {r,[ k |0)) , (5) 

i.e., the eigen state ry lk |0) is also the eigen state of the 
whole graph H . Then such a state represents the trap- 
ping or bound state of a particle within the subgraph I 
with infinite life time. This rigorous conclusion has im- 
portant implications in the design of quantum network 
to store particles in the target region at will. Figure Q] 
represents an arbitrary graph of tight-binding network 
within part of which particles can be confined without 
any leakage. The whole graph can be decomposed into 
two subgraphs I and V which are connected via the cou- 
plings between the joint sites (i,j,h) and (i',j',h'). The 
perfect bound state can be formed in subgraph I as the 
eigen function of Hi when it has wave nodes on all the 
joint sites h). The existence of additional wave nodes 
indicates the multiple bound states can be formed. 
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FIG. 2: (Color online) Configuration of 7r-shaped lattice that 
admits the formation of trapping particles. It consists of three 
chains A, B and C. Two three-site chains A and B are side 
coupled to the chain C with arbitrary number of sites (a). 
This graph can be decomposed into three subgraphs, one of 
which is eleven-site chain enclosed by the dashed rectangular. 
The single-particle eigen states of the eleven-site chain with 
quasi momenta k — n/4,ir/3, and n/2 have two, three and 
four wave nodes, denoted by empty circles in (b), (c) and 
(d). For states with k — n/4 and n/2, the joint sites being 
all wave nodes, the particle can be trapped in the eleven-site 
chain, while state with k = n/3 is not a bound state. States 
with other values of k can be analyzed accordingly. 



III. DEMONSTRATION CONFIGURATIONS 

Now we investigate a class of practical examples to 
demonstrate the application of the result above. We con- 
sider a system of 7r-shaped lattice (Figure ^ , consisting 
of an infinite chain side coupling to two finite chains of 
length No at the joint sites 1 and L, which has the Hamil- 
tonian 

H = H a + Ht, + H c + .Hjojntj (6) 

N 

H a +H b = -K ^(aja i+ i +b\b i+ i +H.c), 

i=l 
oo 

H c = -k (c] ci+i + H.c), 
-Hjoint = -K {a\ci+b\c L + Y[.c.), 

where aj (&j and ct) is the boson or fermion creation 
operator at the jth site in the chain a (b and c). The side 
coupling model was employed to depict coupled-cavity 
system for stopping and storing light coherently [HI . For 
the simple case with the shortest side chains, i.e. Nq = 1, 
the configuration is equivalent to the atom-cavity system 
with single excitation [l(| [H| , where the side-site state 
represents the excited state of the two-level atom. 

First of all, we consider a simplest case: the hop- 
ping integrals are identical for all chains, i.e., k = kq. 
This graph can be decomposed into three subgraphs: left 
chain, right chain and central chain of A = 2Nq + L sites. 
The eigen wave functions of the central chain are given 
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by 
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A + l 



sinfcj, j £ [1, A], 



(7) 



where k — nn j (A + l), n £ [1,A], with corresponding 
eigen values — 2Kcosfc. These states possess the wave 
nodes at 
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(A+l)m 



(8) 



where m are certain integers which ensure the existence 
of integer jk for a given n. Then in the case of that {jk} 
cover the joint sites Nq + 1 and Nq + L simultaneously, 
the corresponding eigen states are the trapping states, 
i.e., a particle can be hold along the central chain for- 
ever. An example of No = 3, L = 5 is depicted in Figure 
where only typical cases with k = tt/4, 7t/3, and w/2 
are presented. Actually, [Eq. shows that states with 
k = 7r/4, 7r/2, and 37r/4 have wave nodes at the joint 
sites. Therefore there are three resonant bound states 
for this configuration. It has been proposed that such 
kind of trapping state can act as a cavity when a boson 
system is considered [ll|. Remarkably, two peculiar fea- 
tures are identified. First, the bound state has infinite 
life time in the ideal case without decoherence since it is 
based on the mechanism of Fano interference rather than 
two potential barriers. Second, the number of the cavity 
mode does not solely depend on the size of the cavity L 
like the case of using infinite potential well for particle 
trapping. For example, taking No = 1, one can achieve a 
single mode cavity with k — ir/2 for arbitrary odd L, but 
none for even L. Meanwhile, it will be shown later that 
there is another type of bound state, evanescent bound 
state. Besides these exact bound states, there exist eigen 
states of the subgraph which have nonzero, but very small 
probability at the joint sites in the case of large L. Such 
kind of state has finite but long life times, which is called 
quasi resonant bound states. To demonstrate these con- 
cepts, we present a numerical simulation of the damping 
process for various modes in two typical systems with 
iVo = 2, L = 4 and N — 3, L — 123, respectively. A 
particle is initially located in the subgraph in the eigen 
states \k) [Eq. (0)]. We investigate the dynamics of the 
states by computing the quantity 
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P(M) = ( !>& + &&)+£ 



c\Ci 



(9) 



i=l 



k.l 



where {■■■) k t denotes the expectation value of the proba- 
bility of the particle within the subgraph for an evolved 
state exp (— iHt) \k). Figure [3] shows the numerical sim- 
ulation of P{k,t) as functions of the mode k and time t 
for a short L in the upper plot while for a longer L in the 
lower plot. There are three types of curves in the two 
plots: (i) remaining unitary; (ii) damping slowly; (iii) 
dropping drastically and then keeping at a finite value. 
Cases (i) occurs in both two configurations, correspond- 
ing to perfect resonant bound states. Case (ii) occurs in 
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FIG. 3: (Color online) The probability P(k, t) of a particle, 
initially located in the state \k), remaining in the subgraph. 
The simulations are performed in two typical systems with 
No — 2, L = 4 (upper panel) and No = 3, L = 123 (lower 
panel) , respectively. The shapes of all the curves can be clas- 
sified into three types as mentioned in the text. 



large- L system, corresponding to quasi resonant bound 
state (we omit such kind of curve in the lower panel). 
Case (iii) occurs in small- L system, corresponding to an- 
other type of bound state, evanescent bound state, which 
will be discussed in detail later. 

A resonant bound-state configuration can be under- 
stood from the point of view of interference. The bound 
state we constructed in this manner is the standing-wave 
like state in the subgraph. In general, the formation of a 
standing wave in a quantum system is due to the infinite 
potential barriers which reflect the wave with any mo- 
mentum. Then there is no additional condition for the 
distance between two barriers. In a tight-binding net- 
work, a side coupled chain can act as the infinite potential 
barriers for the incident wave with certain momentum. 
As an example, it can be readily shown by the method 
below that, for an incident wave with k = tt/2, the trans- 
mission coefficient T through one side coupled chain of 
length iVo can be expressed as T = [l + (-1)^°] /2. It 
has been discussed in Ref. [T3, EH for the case of Nq = 1 . 
Besides the mirror condition T = 0, a matching distance 
between two side coupled chains is also required to form 
a standing wave. This will be discussed below in the aid 
of exact results. 

In the above analysis, the trapping subgraph is the sim- 
plest lattice, an open chain. There are some a little more 
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complicated subgraphs, the hierarchical lattices, as the 
demonstration configurations. It has been shown that 
there are eigen wave functions of these hi- 
erarchical lattices, whose amplitudes are zero at certain 
sites. When these lattices are embedded in a network by 
linking the nodes only, the trapping states are formed. 
Considering an arbitrary generation Vicsek fractal as an 
example, there is eigen wave function whose amplitude 
is zero at the center of every five-site cell. Then when 
such a lattice is embedded in a network by linking the 
center sites only, the corresponding eigenstate is the trap- 
ping state with respect to the network. Nevertheless, for 
the hierarchical lattice itself, this eigenstate becomes an 
extended state as its size grows up. 



IV. BETHE ANSATZ RESULTS 



For simplicity, only simple cases with kq — k are consid- 
ered to demonstrate and explore the obtained rigorous re- 
sults. The existence of the solution requires (L — 1) m — 
{N + 1) n, where n E [1, L - 2] and m e [1, N Q ]. Obvi- 
ously, the resonant bound states in the above mentioned 
example with Nq = 3 and L — 5 is the simplest case of 
m = n = 1, 2, and 3, corresponding to momenta 7r/4, 
7r/2, and 37r/4, respectively. 

For the evanescent bound state, which possesses 
nonzero particle probability at and around the joint 
points, the momenta k and q are determined by equa- 
tions 



-ik(L-l) 



C(A(L-1)) 
E = -2k V (?) 



± 1 



C(?(JVo + l)) 



(14) 
(15) 



We now turn to discuss the complete bound states in a 
subgraph by taking the network of [Eq. ([6])] as an exam- 
ple. It is worthy to point out that the bound states con- 
structed by the before-mentioned method are not com- 
plete. In the following it will be shown that there are two 
types of bound states: resonant and evanescent. The 
former describes trapped particle in a specified spatial 
region and the later describes particle with an exponen- 
tially decaying probability beyond a specified spatial re- 
gion. In the following, we investigate this problem based 
on the Bethe Ansatz approach. Actually, the bound-state 
wave functions of the Hamiltonian [Eq. ©J can be 
expressed as a piecewise function over all sites 



Mj) = 

Mj) 
Mj) 



Cie -iku-i) 

C 2 e ikj + C 3 e- tkl 

C Ae ik{j-L) 

A ie iqj + A 2 e~ iqj 
B ie iqj + B 2 e~ lqi 



for j < 1, 
for 2 -< j -< L, 
for j > L, 

for 1 < j < N , 
for 1 < j < N Q . 



Here ip a ,b,c denote wave functions along chains a, b, and c, 
respectively. The coefficients and momenta 6*1,2,3,4, -A1.2, 
Si. 2, k, and q are determined by matching conditions and 
the corresponding Schrodinger equations [l7| 

^(i + o+) = fi>(j + o-), (10) 

-K j+1J i> (j + 1) - Kj-lJIp (j - 1) = E^j (j) . (11) 

where E is eigen energy, Kj±\j are the corresponding 
hopping integrals. The solutions can be classified in two 
categories: resonant and evanescent ones, which corre- 
spond to zero and nonzero Ci,4, respectively. 

For the resonant bound states, zero Ci,4 lead to zero 
particle probability at the joint points, which is consis- 
tent with the above mentioned rigorous results. In addi- 
tion, the momenta k and q are determined by equations 

sin [fc (£ - 1)] = sin [g (JV + 1)] = 0, (12) 
E = — 2kq cosq = — 2k cos k. (13) 



where C (0) = (e ie - e~ l9 ) /2 and 77 (0) = (e w + e^) /2. 
Taking kq = n, Nq = 3, and L = 5 as an exam- 
ple, we have q = k = 0.382z, n + 0.382z, or 0.19H, 
7r + 0.19H, which correspond to symmetric and antisym- 
metric evanescent bound eigen functions, respectively. 
Furthermore, for the case of kq — k, Nq — 2, and L = 4, 
plotted in Figure[3l we have q = k = 0.382i or 7r + 0.382i. 
Accordingly, three initial states with momenta 7r/9, 27r/9, 
and 47r/9, as well as their counterparts have nonzero over- 
laps with the two evanescent bound states. We are then 
able to obtain the long-time behavior of P(t) as 0.5032, 
0.0027, and 0.0058, which are in agreement with the plots 
in the upper panel of Figure [31 



V. SCATTERING PROBLEMS 

In general, trapping and scattering are two contrary 
phenomena which always refer to localized and extended 
states. In the context of this paper, the resonant bound 
state is essentially standing wave like, consisting of two 
constituents: incident and reflected waves. On the other 
hand, the rigorous result for such bound states has no 
restriction to the size and geometry of the subgraph and 
is applicable to the scattering problem. This is another 
main issue we want to stress in this paper. 

For scattering problem, the input, output waveguides 
and the center system should be involved. One can take 
the input waveguide, which is usually semi-infinite chain, 
together with a part of the center system as the subgraph. 
The resonant bound state in such a subgraph corresponds 
to a total reflection. Actually, the trapping wave function 
within the input waveguide region is the superposition of 
two opposite travelling plane waves with the identical 
amplitudes. They correspond to the incident and total 
reflected waves. And the eigen energy E of this trapping 
state is exactly the transmission zero, i.e., T(E) = 0. 
Taking the above 7r-shaped lattice as an illustrated ex- 
ample, the subgraph containing the input waveguide is 
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depicted by the Hamiltonian 



Wo 



Hi, 



which is a uniform semi-infinite chain. The resonant 
bound states must have a node at site c\ with energy 
E = —2k cos q, where q is determined by the position of 
the node, sin [q (N + 1)] = 0. 

Now we consider the scattering problem of the tt- 
shaped lattice, demonstrating the relation linking the 
scattering state and resonant bound state in the frame- 
work of the Bethe Ansatz. It is worth to note that many 
efforts have been devoted to discuss critically the effect of 
a dangling side coupled chain on the spectrum and trans- 
mission properties of a linear chain, inclu ding the Fano 
resonance, by approximate approaches fill. [l9l I20I [2l| . 

In a 7r-shaped lattice, the scattering wave function has 
the form 



Mj) 
Mj) 



te «(j-i) 



re -ife(i-i) 

+ Be -ik(]-l) 



C b e lq3 + D b e- iq \ 



for j < 1, 
for 2 -< j -< L, 
for j > L, 

for 1 < j < N , 
for 1 < j < N . 



where r and t arc reflection and transmission amplitudes 
for an incident wave with momentum k. Similarly, ap- 
plying the matching conditions [Eq. (fT0|) ] and the cor- 
responding Schrodinger equations [Eq. (jlip j. we then 
obtain 



t 



2 • 2 j 

a sin k 



fc-ia/3sinfc + (/3/2) [ e i2k (L-i) _ i] 



or sm' 



(17) 



where a = Ksin [q (No + 1)] and (3 = Ko sin (qNo). Note 
that zero a leads to vanishing of t, while zero /3 leads to 
vanishing of r. The former and latter are in agreement 
with the conclusions of the above analysis from the inter- 
ference point of view for the total reflection and resonant 
transmission, respectively. 

From Eq. (|17[) . the transmission probability has the 
form 



T 



4 • 4 j 

a sm k 



k + (/3/2) (/? 2 + 4a 2 sin 2 k) sin 2 [k (L - 1) - 5} 



or sm 



(18) 

where tan S — 2asinfc//3. Eq. (fT8|) allows the analytical 
investigation on the transmission features. First, it is 
found that the total reflection condition coincides with 
the resonant bound condition Eq. (fT2|) . It indicates the 
conclusion that an incident wave is totally reflected by 
the side coupled chains if its energy exactly equals to the 
resonant bound state energy. The same conclusion has 
been obtained for some similar systems 0, [ljl H(J Hl[ . 
This is a direct result from the fact that the scattering 
of any dangling side coupled chain is isotropic for the 
incident waves from both sides along the waveguide. On 
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FIG. 4: The plots of T(E) from Eq. JTBJ for the systems of 
iV = 2 (left), 3 (right) with different L. 



the other hand, the resonant transmission condition is 
also easy to be understood from the aspect of wave nodes 
in the subgraph. In fact equation sin (qNo) = indicates 
the effective disconnection of the wave guide from the 
side coupled system. 

Second, for a fixed No, the common transmission ze- 
ros and reflection zeros for arbitrary L can be sim- 
ply determined by a = and f3 = 0, respectively. 
More precisely, for the incident waves with fc m j n 
cos 1 {(k /k) cos \rm/ (N + 1)]}, n € Z, we have T = 0, 
while the one with fc max = cos -1 [(kq/k) cos (nir/No)}, we 
have T — 1. The rest reflection zeros are L-dependent 
and determined by sin 2 [k (L — 1) — 5] = 0. 

The transmission spectra are plotted for n — k , N = 
2 ,3 and different L in Fig. [4] as illustration. We can see 
that the common transmission zeros occur at E = — 1 
for A*o = 2; E — for ^0 = 3, while the common 

reflection zeros occur at E = for No = 2; E = — 1 for 
./Vn = 3, which are in agreement with the above analysis. 
From the plots, one can find that it does not exhibit 
perfect Fano line shape. Nevertheless, the peak and dips 
profiles are the direct result of interference result from 
subwaves in different paths. Actually, the formations of 
fe m i n and fcmax correspond to the complete destructive 
and constructive interferences. 

Now we focus on the L-dependent reflection zeros. 
Consider a system with fixed L = Lo, the L-dependent 
reflection zeros occur at ko, which satisfies 



sin 2 [ko (L - 1) - S\ = 0. 



(19) 



Meanwhile, for a system with L = Lo + in, the corre- 
sponding transmission coefficient obeys 

T(ko,L + m) =T(ko,L -m) (20) 

for Lo — m y 0, due to the identity 

sin 2 [ko (L + m - 1) - 6] = sin 2 (mk Q ) . (21) 

This fact leads to an interesting conclusion. For a cer- 
tain ko, if there are two systems L and L' that sat- 
isfy T(ko,L) = T(ko,L') = 1, there should exist a 
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E (units of k) 

FIG. 5: (Color online) Transmission probability T(E) for the 
configurations with No = 2 (up panel), 5 (down panel) and 
L — 5 (solid blue line), L — 6 (dashed red line). The plots 
show the evident swapping of peak-dip profiles. 

series of different L,L',L",L"',..., satisfy T(k ,L") = 
T(ko,L'") = ... = 1. Especially, applying this conclu- 
sion for m = 1 case, it follows that there is no fco to 
satisfy T (fco, L) — T (fco, L + 1) = 1, except the common 
reflection zeros. In other words, there is no L-dependent 
reflection zeros for L and L + 1 meeting at the same fc. 
This feature enhance the probability of the occurrence of 
so called peak-dip swapping as L changes [2(| Hl[ . 



For a fixed Nq, one can always find two systems with 
successive L, that they have at least one peak (reflection 
zero) located at each side of a common dip (transmis- 
sion zero). Since there is only one peak at each fco, the 
peak-dip swapping profile is formed in the vicinity of a 
common dip. Here we exemplify this point by investigat- 
ing the cases with k = kq and small Nq. For Nq = 2, 
one of the common transmission zero is fc m i n = 7t/3, 
while the L-dependent reflection zeros are determined by 
sin 2 [fc (L — I ) — S] = 0. The closest (or closer) solution 
of fco around fc m i n = 7r/3 are the left one fcoi, = 0.29-7T 
(E = —1.21) for L = 5, and the right one fcofl — 0.36-7T 
(E = —0.84) for L = 6. The profiles of the corresponding 
transmission spectra are plotted in Fig. [5] (upper), which 
exhibit the same character as the one in Fig. 7 of Ref. 
[2fj| . Another example for Nq = 5 and L = 5, 6, is also 
plotted in Fig. [5] (lower). One can see the occurrence of 
the profile of evident peak-dip swapping. 



VI. SUMMARY 

In summary, we show in this paper, within the context 
of a tight-binding model, that a particle can be trapped 
in a nontrivial subgraph. As an application, we examine 
concrete networks consisting of a 7r-shaped lattice. Exact 
solutions for such types of configurations are obtained to 
demonstrate and supplement the rigorous results. It is 
shown that there are two types of bound states: resonant 
and evanescent. We also link the trapping rigorous re- 
sult to the scattering problem for such a subgraph being 
embedded in a one-dimensional chain as the waveguide. 
It is shown that an incident wave experiences total re- 
flection under certain condition. Finally, we also investi- 
gate the scattering features of the 7r-shaped lattice in the 
framework of the Bethe Ansatz. Such rigorous results 
are expected to be necessary and insightful for quantum 
control and engineering. 

We acknowledge the support of the CNSF (Grants No. 
10874091 and No. 2006CB921205). 
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